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Numerical studies together with asymptotic and spectral analysis establish regimes where soliton 
pairs in degenerate optical parametric oscillators fuse, repel, or form bound states. A novel bound 
state stabilized by coupled internal oscillations is predicted. 
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Many nonlinear media can support soliton-likc struc- 
tures when contained in a driven optical cavity [0-0 . We 
will refer to such structures as cavity solitons (CS). In 
quadratic nonlinear media CS have recently been pre- 
dicted in both optical parametric oscillator (OPO) ||] 
and second harmonic generation [H configurations. Al- 
though experimental observation of x CS remains a 
challenge, impressive bistability results |q] demonstrate 
the required level of nonlinearity and thus pave the way 
towards this goal. 

The large values of effective x'-^-' accessible in artifi- 
cially phase-matched materials in combination with their 
practically instantaneous response are important advan- 
tages of using quadratic nonlinearity for implementation 
of CS for all-optical processing of information. They thus 
represent an interesting alternative to the CS which can 
be created in cavities with dispersive-absorptive [E|,and 
resonant electron-hole types of nonlinearities. In all 
such schemes high CS density is desirable and therefore 
understanding of their interaction is a practically impor- 
tant problem which is still largely open. In this Letter 
we focus on the interactions of CS found in the below- 
threshold regime of a degenerate doubly resonant OPO, 
under conditions where the signal field has three coexis- 
tent plane- wave states [y[ . 

Assuming phase-matching, a plane-wave input field, 
and ignoring walk-off, the mean-field OPO equations can 
be presented in the following dimensionless form [0 

-idtEi = [aidl + 5i+ i-ii)Ei + {E2 + ii)El, (1) 

~idtE2 = {a2dl + 62 + il2)E2 + El/2, 

Here Ei and (£'2 + f) are the signal and pump fields, 
respectively, at frequencies w and 2w (we use /.t as a mea- 
sure of the pump strength). The slow time t is scaled so 
that 7m (proportional to the cavity damping rates) and 
5,n (to the detunings from its resonances) are of order 
unity. Here and below m = 1,2. 

This system can describe either diffractivc or dispersive 
effects. We consider x a dimensionless transverse coor- 
dinate, and so set am = ^/m. For this case, existence 
of CS for (5m < was numerically demonstrated 0] for 
^L < A^ < /ifl, w here /il = \^i52 + 72^1 |/\/^i + ll-, and 
Mfl — V^i+7i ^^ *^^ OPO threshold. Within this range 
two different non-trivial homogeneous solutions {E^ ^ 0, 
dxEm = 0) coexist with the trivial one (£',„ — 0), and the 
CS are sech-like localized states on the zero background. 

We start our analysis by applying a perturbative 



method ||6| to the problem of CS interaction. We seek 
solutions of Eqs. (1) in the form 

Em{x, t) = Am{x - Xa) + Bm{x - X b) + 

e{am{x~XA,XB,t) +bm{x-XB,XA,t)) +0{e^), (2) 

where Am{x — xa) and B„i{x — xb) are CS centred on 
xa,b- Note that Eqs. (^ are invariant with respect to a 
■K phase flip of the signal field, so that A and B can be ei- 
ther in-phase or out-of-phase CS. We assume < e <^ 1, 
and that the perturbation functions Om, bm are negligible 
except close to xa, xb respectively. We further assume 
that XA,B vary on the slow time scale t = et and that 
d ^ \xa ~ xb\ is large enough that the overlap functions 
Xi — A2BI + B2AI and I2 = ^iBi are of order e. 

Substituting ansatz ^ into Eqs. (y) and truncating 
O(e^) terms we obtain two analogous systems of equa- 
tions for am and bm, the former expressible in the form: 



{Ca - dt)a = -{drXA)^A +I/e, 



(3) 



Here a ~ (Reai, Rea2, Imai, 777102) ; operator Ca is the 
linearization of Eqs. (13) around the soliton Am] ^a = 
dx{ReAi, ReA2, ImAi,ImA2)^ is the neutral eigenmode 
of £^ associated with translational symmetry, CaCa — 0; 
and I = {—ImXi, —Im'l2, Reli, Rel2Y' controls the in- 
teraction of the two CS. 

The solution of Eq. (||) should in general be expressed 
as a superposition of the eigenmodes ^„ of Ca, C-A^n = 
Ki^n, with time dependent coefficients, because the CS 
interaction will couple to them all. However, apart from 
the above-mentioned neutral eigenmodes, the only ana- 
lytic knowledge about the eigensystems of Ca and Cb 
is that they have two bands of continuum modes with 
eigenvalues A lying on ReX — —"fm, i-e- that all extended 
eigenmodes are damped. We have obtained their full 
eigensystems numerically, using finite-differences, over 
wide ranges of all relevant parameters. We find that for 
sufficiently large dissipation all cavity solitons are stable 
throughout the entire region of their existence. A Hopf 
bifurcation can occur as photon lifetime is increased, but 
we will not consider here any parameter regions where 
isolated CS are unstable. With oscillatory eigenmodes 
absent or well damped, only the neutral mode is easily 
excited by external perturbations, and so we meantime 
neglect all other modes. This enables us to obtain semi- 
analytic results on CS interactions. 



To exclude secularly growing solutions the right-hand 
side of Eq. (g) must be orthogonal to the neutral eigen- 
mode of £^ (which we calculated numerically). This 
solvability condition, together with that for the B soli- 
ton, defines a function / which governs the dynamical 
evolution of the distance d between the soliton centers: 



dtd^fid). 



(4) 



We computed / for both in-phase and out-of-phase inter- 
acting CS, for many parameter values. Typical examples 
are plotted in Fig. Il|. Regions where / is negative (pos- 
itive) correspond to CS attraction (repulsion). Zeros of 
f{d) thus mark stationary bound states of CS pairs, which 
are stable if ddf < where / = 0. 

We find that this equation gives generally correct pre- 
dictions of the inter-soliton forces, in particular that in- 
phase CS attract and out-of-phase CS repel. Both re- 
pulsion and attraction become stronger as fi increases, 
presumably because the signal component {Ei) of the 
CS becomes less localized as /x approaches plane-wave 
threshold at /ifl. A similar effect can be envisaged in 
other CS models. For in-phase CS the function / can de- 
velop pairs of zeros, see Fig. yj. This predicts birth of new 
pairs of CS bound states, one stable and one unstable. 

In Fig. 2 we present simulation results showing differ- 
ent interaction scenarios for two CS initially separated by 
about three soliton widths. First, we consider interaction 
of in-phase solitons. For small fi mutual attraction results 
in fusion of two solitons into one (Fig. 2(a)). Gradually 
increasing /i we first observe formation of a stable oscilla- 
tory bound state (Fig. 2(b)), then of a stationary bound 
state (Fig. 2(c)) which is stable (the radiation visible 
in Fig. 2(c) decays, albeit slowly). Note that the equi- 
librium separation in Fig. 2(c) is predicted quite well 
by the appropriate zero of f{d) in Fig. 1, even though 
these CS are close enough to endanger the assumptions 
of our perturbation method. Stationary two- hump soli- 
tary states have been found previously as solutions of 
an approximate equation derived from Eqs. (El), but no 
analysis of soliton interactions was performed. Note that 
we have found not only two-hump but also multi-hump 
solitary states. However the latter were usually dynam- 
ically unstable. Further details on this issue will be re- 
ported elsewhere. Close to the upper boundary of CS 
existence the interaction of two solitons excites a global 
pattern, see Fig. 2(d), via generation of a switching wave 
from the stable trivial solution up into the modulation- 
ally unstable nontrivial homogeneous state. As predicted 
by Eq. (}^, out-of-phase CS repel each other throughout 
the entire region of their existence - contrast Fig. 2(e) 
with Fig. 2(d), which corresponds to the same value of 

Now we will describe numerical results of the inter- 
action of CS where weakly-damped oscillatory modes 
strongly influence the soliton interactions. Oscillat- 
ing solitons generally radiate energy, which can become 
trapped between neighbouring solitons, exerting a radia- 



tion force which may lead to formation of a bound state. 
An effect of this kind has been reported for solitons in 
models with a weak global dissipation 0. We investi- 
gated a quite different situation, where linear waves es- 
caping from the soliton are strongly damped. Here strong 
interaction between the solitons is due, not to radiation 
modes, but to proto-Hopf modes, and thus has novel as- 
pects. 

The effect is strong providing that two conditions are 
satisfied. First, and crucially, the corresponding eigen- 
modes must have tails with well pronounced and weakly 
decaying oscillatory structure, see Fig. 3(a). Second, 
as might be expected, the oscillatory mode should be 
weakly damped (see Fig. 3(b)), i.e. the CS is close to 
a Hopf instability. If both conditions hold, then, even if 
the global damping due to the 7™ is strong, a CS acts as 
a guide for waves weakly damped in both space and time. 
If a second CS is close enough, these guided waves can 
couple and reinforce each other. Fig. 3(c, d) illustrates 
the dynamics of two interacting CS having eigenmodes 
shown in Fig. 3(a). Note that the separation of the in- 
teracting solitons in Fig. 3(c) is much greater than their 
width. Comparison with Fig. 3(b) clearly indicates that 
the undamped pulsations shown in Fig. 3(d) originate 
from coupling and mutual reinforcement of the oscilla- 
tory modes of the two solitons. A further interesting 
point is that we find these dynamic bound states also for 
out-of-phase solitons, balancing the repulsion induced by 
their neutral mode interaction. 

Quadratic nonlinearity is also known to support soli- 
tons in free propagation geometry and in particularly in- 
teraction of these solitons has recently been studied both 
experimentally and theoretically P,p|. Hamiltonian na- 
ture of free propagating solitons results in their interac- 
tion obeying the laws of Newtonian dynamics 0- ^^~ 
other important difference is that the relative phase of 
the interacting solitons can take only two discrete values 
in a cavity while it is a continuous free parameter in a 
propagation scheme. In spite of these differences fusion of 
the in-phase solitons and repulsion of the out-phase are 
common features in both schemes. However, existence 
of the stationary and oscillatory bound states coupled 
either via translational neutral modes or via internal os- 
cillatory modes are novel important features arising due 
to presence of the external pump and cavity losses. Bal- 
ance between the pump and losses acts as an additional, 
equally important with the balance between diffraction 
and nonlinearity, mechanism of soliton formation inside 
an optical cavity. 

In summary, we have presented the analytical and nu- 
merical study of the interaction of cavity solitons in a de- 
generate OPO and identified distinct static and dynamic 
binding mechanisms. 
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FIG. 1. Plots of the CS velocity function / vs d. Full 
(dashed) lines correspond to in-phase (out-of-phase) solitons 
and thin (thick) lines correspond to /i = 1.6(1.9). Other pa- 
rameters are S2 = —4, Si — —1.8, 71 = 1, 72 = 0.8. 
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FIG. 2. Interaction dynamics of x cavity solitons [9]. 
At different values of pump parameter /i, in-phase CS: (a) 
merge, jj, = 1.6; (b) form oscillatory bound state, /i = 1.8; (c) 
form stable stationary bound state, /i = 1.9; (d) generate a 
pattern via a switching wave, /i = 2. Out-of-phase solitons 
repel, e.g. at jj, — 2, (e). Other parameters as for Fig. 1. 
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FIG. 3. Dynamic interaction of CS for /i = 2, Si = —3, 
S2 ~ —12, 71 = 0.3, 72 = 1, for which the CS has a mode 
with eigenvalue pair A ~ —0.03 ± i4.14: (a) Spatial structure 
of the eigenmode, Re{ui) - full lines, Re{u2) -dashed lines; 

(b) Temporal evolution of signal energy Q — ^ dx\Ei\^ for 
slightly perturbed single soliton, showing damped oscillation; 

(c) Spatio-temporal evolution [9] of \Ei\ (time window much 
later than in (b)), showing dynamic bound state; (d) Tem- 
poral evolution of signal energies of the two CS in (c) , in the 
same time window, showing rapid undamped oscillations and 
slow energy exchange between the two solitons. 



